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1. INTRODUCTION 

Let if be a differential operator with quadratic form 

Q(f) = I E a i , j (x)D i f(x)WJ^)d 

w I -I j> 



\i\<m 
\j\<™ 



where aij(x) = aj,i(x) are complex- valued bounded measurable functions 
on a region Q C WL . Dirichlet boundary conditions are imposed upon H 
by restricting the domain of the quadratic form to be the Sobolev space 
W™' 2 (tt), which is the closure of C c °°(fi) in the Hilbert space W m > 2 {Vt). Here 
C^°(f2) denotes the space of smooth, compactly supported functions on Vt 
and W m ' 2 (Q) is the space of all functions / G L 2 (Q) whose weak derivatives 
D a f lie in L 2 (fl) for all multi-indices a such that \a\ < m. It is equipped 
with the inner product 

(f,g)m?:= E {D a f\D a g) 2 . (2) 

|o|<m 
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In the special case where the coefficients are constant and chosen so that 

\i\ = \j\=m 

and ciij = whenever |i| + \ j\ < 2m — 1, the associated operator H =: i?n,m 
is the polyharmonic operator (— A) m |oiR of the region Q. We denote this 
particular quadratic form by Q m . 

Throughout this paper we shall assume that the coefficients 
chosen in such a way that Q satisfies the Garding inequality 

XQm(f) - 4f\H < Q(f) < VQm(f) + 4f\g (3) 

where < A < fi and c, d are non-negative constants. We then say that H 
is uniformly elliptic. Quadratic forms satisfying the Garding inequality are 
closed on the domain W m ' 2 (Q). For a more detailed account of uniformly 
elliptic operators see J3J. 

Note 1.1. If the sum in equation ([!]) is only taken over non-negative multi- 
indices i, j with \i\ = \j\ = m then the operator is said to be homogeneous 
of order 2m. If this is the case then we may set c = d = in the Garding 
inequality ([$[). 

LEMMA 1.2. For N < 2m the operator H has a heat kernel K(t,x,y) 
which satisfies 

\K(t,x,y)\ (4) 

for all ijeS] and all t > 0. 

Proof. See [f|, Corollary 15 and Lemma 17]. If the operator H is homoge- 
neous then the bound is valid even without the term e*. □ 

In 0, Lemma 19], Davies obtains the pointwise heat kernel bound 

\K(t,X,y)\ < Cl r N ' 2m exp[-C 2 |y - x |2m/(2m-l) f -l/(2m-l) + k ^ ( 5 ) 

for all t > and x,y G ~R N , where ci,C2,k are positive constants. If the 
operator is homogeneous of order 2m then the constant k may be set to zero. 
In 0, Barbatis and Davies find the sharp constant c 2 for this expression in 
the following sense: If H is uniformly elliptic, homogeneous of order 2m, and 
satisfies the Garding inequality 

Qm(f) < QU) < VQm(f) (6) 
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then 



\K(t,x,y)\ < c e r N/2m exp[-(a m - 0(fi - 1) - e)\y - x \^/(2m-i) t ~i/{2m-i)^ 

(7) 

where e > and 

a m = (2m - l)(2m)- 2m/(2m - 1} sin[7r/(4m - 2)]. (8) 

This is shown to be sharp by considering the case H = H R N jrn . 

The Euclidean metric do(x,y) := \y — x\ is a relatively weak and unnat- 
ural way of expressing the heat kernel bound for non-convex regions. As 
an example, for a horse-shoe shaped region whose extremities are touching, 
the Euclidean distance do(x,y) can be made arbitrarily small for internally 
distant points x, y G Q, rendering the heat kernel bound useless at these 
points. 

The aim of this paper is to find heat kernel bounds which are given in terms 
of the Riemannian metric d g (see definition (|12|)), instead of the Euclidean 
metric do- This would improve the original bound; for a horse-shoe shaped 
region which touches itself or even overlaps with itself one may choose x and 
y to make do(x, y) arbitrarily small whilst d g (x, y) remains large. In order to 
find bounds involving the Riemannian metric, we first find bounds involving 
Riemannian-type metrics. 

Definition 1.3. For (3 > define the Riemannian-type metrics d m $ : 
n 2 -> R + by 

d m ,p{x, y) = sup{0(y) - <f>(x) : <p e £ m ,p\, (9) 

where denotes the set of all bounded real valued smooth functions (ft on 
a region Q such that 

||V0|U<1 and IpVlloo < (10) 

for all non-negative multi-indices i such that 2 < |i| < m. 

We show, for an arbitrary uniformly elliptic operator whose coefficients 
need only be measurable, that 

\K(t,x,y)\ < c 1 r Ar / 2m exp[-c 2 d m/3 (x,?/) 2m/(2m " 1) r 1/(2m - 1) + k{l + p 2m )t] 

(11) 

where iV < 2m. Note that since d m< o is the Euclidean metric do, setting 
(3 = in this equation retrieves equation (|5|). 
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If Q is a region with C 2 boundary and radii of curvature uniformly bounded 
below by r then there exists a constant K, dependent only upon m and N, 
such that for (3 > 4K/r, 

where d g is the standard Riemannian distance 

d g (x,y) := inf{Z(7) : 7(0) = 7(1) = y, 7 C f2, 7 is cts and piecewise C 1 }, 

(12) 



and where ^(7) is the length of the path 7. See Theorem 4.18 . This feature 



allows us to convert the off-diagonal bound ([□]) into the bound 

\K(t,x,y)\ < c 1 r N/2m exp[-c 2 d g (x,y) 2m/{2m - 1) t- 1/(2m - 1) + kt] 

(13) 

involving the Riemannian metric. 

Note 1.4. As a special case of [Rl Theorem 3.2.7 and Corollary 3.2.8] 
Davies has the bound 

< K(t, x, y) < c 5 t~ N / 2 exp{-d g (x, y) 2 /A(l + 5)t}, (14) 

valid for all N, on the heat kernel K(t, x, y) of the Dirichlet Laplacian on a 
region Q C Mr. The regularity of the boundary of Q, required for higher 
order operators is a genuine feature of the order. □ 

A more natural class of metrics for the determination of heat kernel bounds 
is, in a certain sense, the class of Finsler-type metrics d a ^ M induced by the 
operator itself. Here a denotes the principal symbol 

a(x,£)= £ a itj {x)C +j (15) 

\i\=m 
\j\=m 

of the operator, and M is a positive constant. 
Definition 1.5. The Finsler-type metrics are defined by 

d a ,M{x,y) = sup{0(y) - 4>(x) : G #o,m}, (16) 

where JP a ,M denotes the set of all bounded real-valued smooth functions on 
Q such that 

a(x,V0(x))<l and H-D^Hoo < M (17) 
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for all non-negative multi-indices i such that 2 < |i| < m. 

Note that the Riemannian-type metrics d m ^ are similar to the Finsler-type 
metrics induced by the polyharmonic operator (— A) m |Dm, whose symbol is 

a{ X) = \i\ 2m - 

Under the following assumptions, Barbatis |I| uses Finsler-type metrics to 
express a sharp heat kernel bound. 

Assumptions 1.6. 

(i) H is uniformly elliptic and homogeneous of order 2m > N; 

(ii) The coefficients aij lie in the Sobolev space W m '°°(fl); 

(iii) The symbol of H is strongly convex (see Definition |1.7| below). 

Definition 1.7. For \k\ = 2m define 

A;' 

^ '' H=lil= OT 

i+j=k 



We may rewrite the principal symbol of H as 

Im 



a(x,0 = { ^T^ a k{x)i k 



\k\=2m 

We say that the symbol a(x, £) is strongly convex if the quadratic form 

r(^, C) = z2 a p+q( x )CpCq 

\p\=m 
\q\=m 

is non-negative for each x G fl, where ( = (C P )| P |=m £ ^ and 

fn + m — 1\ 
V= V n-1 ) 
is the number of distinct multi-indices p with \p\ = m. 

Barbatis proves that, under these assumptions, we have the heat kernel 
bound 

\K(t,x,y)\ < c s t- N / 2m exp[~(a m - 5)d aM ( x , y)^/^)^ 2 ^} 

(18) 

for M large and for t/d a> M < Tg t M- We shall apply this result under the same 
assumptions, to find the sharp constant c 2 in the bound (P~5|). 
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2. THE HEAT KERNEL BOUNDS 



We shall use the technique found in [|J of twisting the operator H to define 
H a 4> = e a ^He~ a( ^, which has quadratic form 



QM) = j n E a,„ J (x){e^ J D i e-^/(x)}{e-^e^/(a;)}d 



N x. 



W<™ (19) 

\j\<m 

Here a > and G <^ m ,/3 for some (3 > 0. The following proposition is a 
generalisation of j|, Lemmas 1 and 2], treating derivatives of the function 
more delicately. 

Proposition 2.1. The twisted quadratic form satisfies the inequality 

\Qa*(J) - W)| < eQ(J) + Ce(l + ^ + P 2m )\\f\\l (20) 

where e may be taken arbitrarily small. 
Proof. Each term in Q a cj>{f) may be expanded using formulae of the type 

e a<t> D i e -*4>f = £>*/ + £ c *( II aD k ^)D ko f, (21) 

r=l 

where the sum is taken over all integers p and non-negative multi-indices 
k , . . . ,k p such that 

p 

'^ / k r = i, k ^ i, and ki,...k p ^0. (22) 
Combining these terms, 
Q*<l>(f)=Q(f)+ I Y,' c kA^ P (f[D k ^)D ko f.^(f[D l ^)D^d N x, 

n r=l «=1 (23) 

where the sum J2' is taken over all integers p, q and non-negative multi-indices 
k , . . . , k p , l , . . . ,l q such that 

p i 

2_, \kr\ < m , E \h\ < m > an d I Ato I + \h\ < 2m — 1. (24) 

r=0 s=0 

Hence 

- q(/)i < e' iicMiiooii« p ( n D^)D*»/ii 2 iia«( n ^0)^7112. 

r=l s=l (25) 



Now 

p 



\a p 



(n^)^°/ii2 < (n ii^iioc)n^ fco /ii2 

r=l 

< (f[(3^- l )\\a p D k0 f\\ 2 



r=l r=l 

V 



r=l 

|| a P/3»-l'-|*ol£)*o/|| 2) 



where s := X)r=o K < m, and so 



\a p > 



r=l 



Jm. n 

< J rn [e|£| 2 * + c e (a 2s + I3 2s )} | 2 d^ (26) 

< / [e|e| 2m + c e ( a 2 - + /3 2w + l)]|/(0| 2 d^ 



(27) 

where e may be arbitrarily small in (26) provided \k \ < s — 1 and arbitrarily 
small in (27) provided \k \ < m — 1. Each term in J2' nas & t least one of 
\ko\ < m — 1 or | l Q \ < m — 1 and so is dominated by 

eQ m (f) + c,(l + a 2m + (3 2m )\\f\\l 
Using the Garding inequality (^) this in turn is dominated by 

eQ(f) + c e (l + a 2m + p 2m )\\f\\l □ 



Note 2.2. If H is homogeneous (see Definition |1 . 1| ) , then 2s = 2m and 
the above proof finishes with equation (26) instead of (27). This yields the 
twisted form inequality 

\Qadf) - W)l < + Ce(« 2m + (3 2m )\\f\\l (28) 

instead of equation Q2"0p. This will induce a corresponding change in inequal- 



ities d|), (H) and @). □ 
Lemma 2.3. There exist positive constants c,k>0 such that 



-H a( f,t 



<exp[A;(l + a 4 + /? 4 )t] (29) 
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and 

WH^e-^W < cr 1 exp[k(l + a 4 + /5 4 )t]. (30) 

for all t > 0, a, (3 > 0, and <p G £ m $. 

Proof. See @, Lemmas 6 and 7]. □ 
Lemma 2.4. Let d m $ be the Riemannian-type metrics of Definition |1.3| . 



There exist positive constants Ci,C2,k such that 

\K(t,x,y)\ < c l t- N / 2m exp[-c 2 d m( ,(x,y) 2m/i2m ' 1 h^ 2m ^ + k(l + [3 2m )t] 

(31) 

for all (3, t > and all x, y G Q. 

Proof. Put ft = exp[— H a ^i\f . Then 

< c ||(-A) m / 2 /t||f /2ni|U!||1_jV/2m 



t||oo ^ HK~ ^ Jt\\2 II J*ll2 

<cg(/,) JV / 4m ||/,nr 7V/2m 

<c{ReQ^(/^ + (l + a 2m + /3 2w )||/ t ||i} JV / 4 -||/ t ||^ iV/2m 

< c{||^/ t || 2 ||/ t || 2 + (1 + a 2m + /? 2m )||/ i |l2} JV/4m H/*ll2~ 7V/2m 

< cjr 1 + (1 + a 2m + (3 2m )} N / Am exp[Jfe(l + a 2m + /3 2m )t] ||/|| 2 

< ct^ 4 " 1 exp[ife(l + a 2m + (3 2m ^ 



Therefore 

|| expf-^IU^ < ct- Af / 4m exp[A;(l + a 2m + (3 2m )t}. (32) 

By duality, 

\\eM-HaMoo,! < ct- N / 2m exp[k(l + a 2m + p 2m )t]. (33) 
But exp[— H a ^t\ has kernel 

K aHx) (t,x,y) = e a ^K(t,x,y)e- Q ^ (34) 
for allt > and x,y G Q. Equivalently, 

\K{t, x,y)\< ct~ N/2m exp[a((f>(y) - <j>(x)) + fc(l + a 2m + /3 2m )t]. 

(35) 

Taking the infimum over all <fi G i n this bound, we see that 

\K(t,x,y)\ < ct~ N/2m exp[-ad m ^{x,y) + k(l + a 2m + (3 2m )t}. 
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Optimising with respect to a gives 

\K(t,x,y)\ < cr JV/2m exp[-A;'d m ^(x, 2 /) 2m/{2m - 1) t- 1/(2m - 1) + k(l + (3 2m )t\. 

□ 



Theorem 2.5. Let Vl c R N be a region whose boundary is C 2 , with radii 
of curvature uniformly bounded below, and suppose that N < 2m. There exist 
positive constants c\, C2, k such that 

\K(t,x,y)\ < cit JV/2m exp[-c 2 d fl (x,?/) 2m/(2m - 1) t- 1/(2m - 1) + kt] 

(36) 

Proof. This follows by setting (3 = AK/r and applying the main result, 
Theorem gig , of Section |. □ 



Note 2.6. The statement ( p6|) is equivalent to the existence of positive 
constants ci, c 2 , T such that for t < Td g (x, y) 

\K(t,x,y)\ < c 1 t' N / 2m eM-C2d 9 (x,y) 2m ^ m - 1 k- 1 ^ m - 1 ^} (37) 



See Lemma 3.1. □ 



3. SHARP CONSTANTS FOR THE HEAT KERNEL BOUND 

Lemma 3.1. Let H be a uniformly elliptic operator acting in L 2 (Q) where 
Q C M. N and N < 2m. Let K(t,x,y) be the heat kernel of H , and let c 2 be 
fixed. The following conditions are equivalent: 

(i) For all e > there exist positive constants c±,T such that 

\K(t,x,y)\ < c 1 t- 7V / 2m exp[-(c 2 - e)d g (x,y) 2m ^ m - 1 h- 1 ^ 2m ^} 

(38) 

for t/d g (x,y) < T. 

(ii) For all e > there exist positive constants Ci, k such that 

\K(t,X,y)\ < Ci r N / 2m exp[-(c 2 - e)d g (x, y)^/(2m-l) t -l/(2m-l) + fcf ] 

(39) 

for t > and x,y e O. 
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Proof. Suppose that (i) holds. Let k = (c 2 - e )T- 2 W(2m-i) + i. Then for 
t/d g > T 

_ ( C2 _ e ) d 2m/(2m-l) r l/(2m-l) > + t 

and hence using Lemma [L~2", 



|^(t,x,3/)| < cT^^expfl 

< Ct-^ 2m exp[-( C2 - e ) d 2m/(2m-l) t -l/(2m-l) + fcf ] 

Conversely, suppose that (ii) holds. Let T = (e/k)^ 2m ~ l ^ 2m . For t/d g (x,y) < 
T we have fct < e rf2m/(2m-i) r i/(2m-i) and SQ 

\K(t,X,y)\ < Cit- JV / 2m eX P [-C 2 4" I /(2— D r l/(2m-l) + ed 2m/(2m-l) t -l/(2m-l) ] 
= C i r Af / 2m eX P [-(c 2 - e ) d 2m/(2m-l) r l/(2m-l)]_ g 



Theorem 3.2. Suppose that Q is a region with C 2 boundary and whose 
radii of curvature are bounded below. Suppose also that H satisfies the as- 
sumptions |1.6|. For e > there exist positive constants c e , k e such that 



\K(t,x,y)\ < c e t-^ 2m exp[-(a m - e^-V^d^y)^'^ 1 ^ 1 '^ + k e t]. 

(40) 

Equivalently, for e > there exist positive constants c e , T e such that 

\K(t,x,y)\ < c £ t- N/2m exp[-(a m - e)^- 1 ^ aro - 1 ^ fl (a:,y) am/(2m - 1) r 1 ^ !!,n -^] 

(41) 

for t/d g (x,y) < T e . 

Proof. For e > fixed let (3 > AK/r be large enough such that the 
result (HD of Barbatis is valid for M := /i" 1 / 2 ™/?™" 1 , and such that 

(a m - e) < (cr m -5)^1- 

The Garding inequality implies that the symbol a(x, £) defined in equa- 
tion ( [l~5|) satisfies 

A|e| 2m <a(x,0</^le| 2m 
for all (GM^. Let G C/3 and define ^ = /i" 1 / 2 ™^ Then 

a(x,V^(x)) =^ 1 a(x, V0(x)) < |V0(x))| 2m < 1 
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and 

PVIU = jr 1/am ||£>Vl|oo < ^ 1 ' 2m (3^~ 1 < M, 

so ^ G ^ a ,M- Hence 

/i- 1/2m ^(x,y) = sup^- 1 / 2 ™^) - /i- 1 / 2 " 1 ^^) : G C/?} 
< sup{?/%) - ip{x) : V e ^ ,Af} 
= d aM (x,y). (42) 



For t/d g (x,y) < (1 — JK/flr)/^ l / 2m T$ } M ='■ T e we see, using equation (42) 
and Theorem [4.18| , that 

t/d a>M (x,y) < fi 1/2m t/d mif3 (x,y) 

<^ 2m {l-jK/P~ryH/d g {x,y) 

and so using (JT~5|) , 

\K{t,x,y)\ 

< c s t N / 2m exp[-(a m - 8)d a M^y) 2m/(2m - 1) t- 1/{2m - 1) ] 

< c s t N ' 2m exp[-(a m - 5) / i- 1 /( 2m - 1 )rf m , /3 (x, 2/ ) 2m /( 2m - 1 )t- 1 /(^-D] 

< c s t N / 2m ex V [-(a m - (J)//- 1 /^- !) (\ - ^\ d g {x )y ) 2m l {2m ' l h- l ^ 2m ^} 

< c e t N / 2m exp[-(a m - ^-VIM^^jW^-ilri/^-i)], 



Using Lemma [O] this is equivalent to the bound ( f40|) . □ 



4. A COMPARISON OF PJEMANNIAN-TYPE METRICS WITH THE 
STANDARD RIEMANNIAN METRIC FOR A HIGHLY NON-CONVEX 
REGION 

The purpose of this section is to prove a geometrical result concerning met- 
rics, which can be applied in Section |^. The Riemannian-type metrics of 



Definition |1.3| are used in Section ^| to express certain Gaussian heat kernel 
estimates. It is beneficial to compare these metrics for different values of (3 
in order to convert the bound into terms of the standard Riemannian metric 
on a region. The comparison is particularly interesting when one notes that 
d m fl = d is the standard Euclidean metric 

do(x,y) = \y - x\, (43) 



11 



and c? TO)00 is the standard Riemannian metric d g , defined in equation (|12"D. If 
Cl is convex then all the above metrics are identical. For non-convex regions 
however this is not the case, and a useful comparison is non-trivial. Clearly 
d m $ is an increasing function of (5. 

Let Cl be a region in M. N , whose boundary dCl is C 2 , with radii of curvature 
uniformly bounded below. We shall prove that for (3 > AK/r, 



1 




(44) 



where r > is the greatest lower bound of the radii of curvature of the 
boundary of Cl, and K depends only on m and N. This result is also valid 



for locally Euclidean Riemannian manifolds. See Note [4.19 



We now develop tools for local representation of the surface dCl of Cl. 

Definition 4.1. We say that a region Cl C has C 2 boundary dCl if 
dCl = dCl and if for each point p e dCl, there exist a set U = U(p) open in 
M N and containing p, a local coordinate system y = . . . , j/jv-i) and y^, 
with (y, y N ) = (0, 0) at p, and a function h = h(.,p) such that dCl R U has a 
representation 



where G 



Vn = h(y) yeG, he C 2 (G), 
and convex. 



(45) 



G(p) is open in 

The surface dCl within U may equivalently be represented by the non- 
degenerate, bijective C 2 map a : G — > dCl PI U defined in local coordinates 
by 



z(y) = (y,*(y)). 



(46) 



Notation 4.2. 

(i) The matrix derivative of a function f : M. N ~ X — > IR M zs defined by 

( dfj_ dh dfM_ \ 



dy 



9^1 
dfi 



9y 1 
df M 



(47) 



w/iere fa := 71* o / (i 



M) are £/ie coordinate functions of f ; 
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Figure 1: Representation of the boundary 

(ii) The second derivative of a function h : M. N ~ 1 — > H. is the matrix defined 
by 



d 2 h d fdh 



/ d 2 h 
I dy* 



d 2 h 



dy 2 ' dy V dy 



dy 1 l 



(48) 



d 2 h d 2 h 
\dy N _ 1 dy 1 ' ' ' dy 2 N _ 1 ) 



Let the y^-axis point into Q. The unit normal n = n{p) to the surface at 
p — <r(y) is defined in local coordinates by 



dh 



T 



dh T dh^ 



The C 1 map r : G x 



defined by 
r(y,u) = <x(y) + un(y) 



(49) 



(50) 



is non-degenerate at (0,0). 
Define vr : dVl x R -> <9fi x R by 

7r(p, m) = (p, 0) p6 <9f2, u G 

and p : dtt x R -> M w by 

p(p, u) = p + un{p) p G dQ, u G 



(51) 



(52) 



Definition 4.3. For 5 > we define the 5 -neighbourhoods of dQ and O 



by 



(dQ) 5 := {z G M. N : dQ) < 5} 



(53) 
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and 



{zeR N : d(z,n) < 5}. 



(54) 



Proposition 4.4. The 5 -neighbourhood (dQ)s of the boundary is the im- 
age of dQ x (—5, 5) under the map p. Similarly Q$ = ^ U p(dfl x [0, 5)). 

Proof. Suppose that z = p{p,u) = p + un(p) for some p G dQ and some 
— 5 < u < 5. Then \z — p\ = \u\ < 5 so d(z, dfl) < 5. 

Conversely, suppose that < d :— d(z,dQ) < 5. Let p G dQ be such 
that \z — p\ = d. Using the representation of the surface cr(y) = (y,h(y)) 
in the local coordinate system based at p we see that the function y i— > 
\z — (y, h(y))\ 2 is minimized at y = 0. Thus for i — 1, . . . , N — 1, 







yz-(y,h(y))\ 2 



0,1,0, 




y=o 



2{&-(y,h(y)),- 
-2 (z- (0,0). | 0. 1.0 



so z — p is normal to the surface. The vector (z — p)/d has unit modulus so 
(z — p)/d = ±n(p), the sign being dependent on whether {z— p)/d is inward 
or outward pointing. Thus z = p ± cfn(p) G p(<9f2 x (—5, 5)). 

The proof that Q$ — ^ U p(9f2 x [0, 5)) is similar. □ 

Condition 4.5. Let Q be a region in M. N with C 2 boundary such that there 
exists an r > whereby 



B(p±rn(p);r) D <9fi = 



(55) 



/or a// p G 

This condition is slightly stronger than requiring that the radii of curvature 
at points of the boundary are bounded below by r. This is done to exclude 
certain regions for which the results of this section still hold, but which 
require a more technical treatment. See Note |4.19| . 

Lemma 4.6. Equation (|55D in Condition is equivalent to bijectivity of 
the restriction p : dfl x (— r, r) — ► (dQ) r . 

Proof. Suppose that for some p G dQ, B(p + rn(p);r) fl dQ ^ 0. By 
writing B(p + rn(p); r) = Uo<d<r Bd where B d := B(p + dn(p); d) we see that 
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Bd PI dVt ^ for some d < r. Let x G -B^ fl cftl Then d(x,p + dn{p)) < d so, 



by Proposition |4.4j , p + dn{p) G (cKT)^. Hence 

p{p,d) G p(<9fi x (-d,d)), 

and p is not injective. The same conclusion is drawn if B(p—rn(p); r)C\dQ ^ 
for some p G <9fi. 

Conversely suppose that the restriction of p is not injective. Then 

p + un{p) = q + vn(q) (56) 

for some 

(p,u) (q,v). (57) 

Without loss of generality we may assume that \v\ < \u\. Assume also that 
< u < r. Then 

\p + rn(p) — q\ = \vn(q) + (r — u)n(p)\ 

< \v\ + r — u 

< r. (58) 

Moreover, using inequality (|56| ) we see that (0) implies that un{p) ^ vn(q) 
and hence the inequality is strict. Thus 

B(p + rn(p);r)ndQjt®. (59) 

If we assume that — r < u < then similarly we obtain 

B(p - rn(p);r) n dQ ^ 0. □ 



Proposition 4.7. IfQ is bounded then it satisfies Condition fL5 . 

Proof. Since r is non-degenerate at (0,0), an application of the inverse 
function theorem shows that r it is injective and non-degenerate in an open 
neighbourhood iV of (0, 0). The set t(N) is open because r is non-degenerate, 
and since r is continuous, r~ 1 (r(A r ) fl U) is open. See Figure |[ If G and 
r' > are small enough, 

G x (-r',r') C t-\t(N) HU) C N. 

Again, since r is non-degenerate in N, the coordinate neighbourhood V := 
t(G x (— r',r')) is open and so there exists r" > such that B(p; ^r'p) C V. 
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Figure 2: t(N) for a typical region Q. 



Figure 3: Construction of 
Gx(-r',r'). 



See Figure [3|. The collection {B(p; \r'p) : p G <9fi} forms an open covering 
of dQ and has a finite subcovering {B{pi, Ir")}™^ by compactness. Let 



.,m I i ■ 



3' i ) 



This construction has been chosen so that if p G <9f2 then p G 
for some z and then -B(p; 2r) C rf ) C V^. 

We shall now prove that p is injective. For suppose otherwise, then p(p, u) = 
p(q, v) where p,q G dQ, —r<u,v<r. By the construction above, B(p; 2r) C 
Vi for some i. Since p + im(p) = q + vn(q) we see that \p — q\ < \u\ + \v\ < 2r 
so g G V;. Thus p, q & Vi P\ dQ and so 

P = Ti(x,0) = CTj(x) 

9 = n(y,o) = o-i(y) 

for some x, y G Gj. Now Tj(x, w) = p(p, u) = p(q, v) = Tj(y, v), and since Tj is 
injective we see that x = y and u = v. Moreover, p = <Xj(x) = crj(y) = q. □ 

From this point onwards we shall assume that all regions satisfy Condi- 



tion 4.5 



Definition 4.8. We say that a real symmetric (N — 1) x (N — 1) matrix 
A is non-negative, and write A > i/cL4a T > for all a G M^ -1 . 

Proposition 4.9. Let p G <9f2 and let h : G — > K. fte £/ie representation of 
the surface dQ in the local coordinate system (y,yN) based about p. Then 



dhdh T d 2 h 

In-x + -r~r ^ rn ^T^' 

dy dy dy^ 



n 



d/id/i T > 
Ar_1 dy dy 



-rnN 



d 2 h 



(60) 
(61) 
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where 



dh T dh^ 1/2 



\ + dy dv 



(62) 



is the N-th entry of the unit normal n(y) in local coordinates. 

Proof. Let d_{y) = (y, <?(y)) represent, in local coordinates, the surface of 
the sphere, radius r touching the surface (y, h(y)) at the point q = (yo, h(yo)) 
where the unit normal is n = (n, n^)- Then 

(£(y)-g)(£(y)-2™- g ) T = o. 



(63) 



Differentiating with respect to y = (yi, . . . , yjv-i) we see that 



O 1 = \ In~i 
= 2 \In-i 



dg_ 
dy 



(a(y) - 2rn - qf + I N ^ 



dy 



(*(y) - q) : 



^ j (cr(y) - rn - ?) T 



2(y T - rn T - y T ) + 2(# - rn N - h(y )) 



dg_ 
dy' 



Differentiating the transpose, 







l N-l 



+ 



dg dg 



+ (g- rn N - h(y )) 



d 2 g 



dy dy ' ' *•'"' "dy 2 

Since the sphere touches the surface at y = y , we see that 

d 2 g 



(64) 



(65) 



g(yo) = h(y ), ^(y ) = ^ 



(yo), and ^j(yo) < 



d 2 h 
dy 1 



(yo) 



Thus, at y = y , 



l N-l 



+ 



dhdh T 

dy dy 



! Af-l 



+ 



dg dg 
dy dy 



d 2 g d 2 h 
m N —^ > rn N - 



dy 2 



dy 2 



(66) 



(67) 



This result holds for all y € G. The proof of part (ii) uses the fact that a 
ball of radius r fits inside the region. □ 

For a G M^ 1 define 



a 



' dh T dh" 
dy dy, 



-V2 , 2 



d 2 /i 



dn dn T T 
a — a 



a I N -i + 



dhdh 
dy dy 



dy dy 



a \ lN -^Tyd-y K 



(68) 
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and let A be the real symmetric matrix 



Lemma 4.10. 

.. . dn dn T 
(i) 



, T dhdh T 
A = I N -i + -T--T ■ 
dy dy 



\ + dhT(ih \ 1 d 2 h A _ 1 d 2 h 



(69) 



dy dy \^ dy dy J dy 2 dy 2 ' 

(ii) a 2 ,13 G [0,1/r 2 ] ; 

(iii) 1 + 2ar + (3r 2 > and 1 - 2ar + /5r 2 > 0; 

\u\\ ' 



(iv) 1 



I u I 

r 



< 1 + 2cm + (3u 2 < ( 1 + — J for all ue [-r, r]. 



Proof, (i) By differentiating the expression for the normal, we see that 



dn 



1 + 


d/i T d/A 


" 3/2 d 2 /i 


dy dy; 


dy 2 


1 + 


d/i T d/A 


- 1/2 d 2 h 


dy dy j 


dy 2 



d/i T d/A / . „ 

I +— — I (In-iIO 1 



dy dy J 
(In-i\O t 



dh (dh l 
dy Vdy 



^ d/i T d/A 1 d/i {dh r N 
dy dy/ dy \dy ' 



Hence 



dn dn 



d/i T d/z\ 1 d 2 /i 



dy dy 



dy dy I dy 2 



! JV-1 



/ dh T dh\ 1 d/id/i T 
V dy dy/ dy dy 



dh T dh\ dh T 



( d/i T d/A 2 dh( 1 __ 
I dy dyi dy V dy dyi dy 



dT 2 



dfr T d/A 1 d 2 h 



-T- 



dy dy/ dy 2 



; _d/id/i T / d/i T d/A 
JV ~ 1 dydy I dy dyi 



-i 



d 2 /i 
dy^ 



dh T dh\ d 2 h A ,d 2 h 
1 + — — -^A' 1 — . 



\ dy dy J dy 2 dy 2 
(ii) Proposition |4J^ immediately implies that 

1 1 

--<«<-. 
r r 
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Since A is a positive definite matrix, we may pre- and post-multiply the 
results of Proposition [O] by A~ x l 2 to get 



7 JV _ 1 + mj vA- 1 /^A- 1 ^>0; 



dy 2 ' 
d 2 h 
dp 



(70) 



- rn^-^A^ > o. (71) 



The matrices in inequalities (70) and (71) commute, so their product is also 
positive. Hence 

iN-.-r^A-^A-^A-^yO. (72) 

By pre- and post-multiplying by A 1 / 2 , we see that 

( T dhdh T \ 2 2 d 2 h . ,d 2 h n 



so using part (i) 



dndn T 2 2 d /i x d /t T 
dy dy ^ dy 2 dy 2 



Thus /3 < 1/r 2 . 

(iii) Squaring the left hand side matrix of inequality (70) will yield a posi- 
tive matrix 

I N ^ + 2rn N A-^ 2 A^ 2 + r 2 n 2 N A-^A-^A^ 2 . 
Pre- and post-multiplying by A 1 / 2 , we see that 



dhdh T \ r-r d 2 h x „ „ d 2 h , , d 2 h ,,. 

a Jjv-i + 3-3- U T + 2r^a— a T + r 2 ^a— A" 1 — a T > 
\ dydy J dy 2 dy 2 dy 2 

for all a e M^ -1 . Thus 1 + 2ar + (3r 2 > 0. Squaring the left hand side of 
inequality (71), we see that 1 — 2ar + (3r 2 > 0. 

(iv) By part (iii), the polynomial 1 + 2au + f3u 2 dominates 

u u 2 
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at u — 0,r. Moreover, since (3 < 1/r 2 , this holds true for all u G [0, r\. Thus 

1 + 2au + f3u 2 > (l - ^ 
for all u G [0, r]. Similarly, 



1 + 2cra + /fa* > ( 1 + - J 



for all u G [— r, 0]. Also, by part (ii), 



u\ \u\ 



1 + 2au + (3u z < 1 + 2— + 



1 + 



□ 



Let P : G x (— r, r) -^Gx (— r, r) be the projection defined by P(y,u) 
(y,0). Let 



r'(.,.):Gx(-r,r)^M N 
denote the Jacobian matrix 



(74) 



dr 



d(y,u) 



of r. 



Proposition 4.11. Let (y, u) e G x (-r,r) and Ze£ n = (a, a N ) e R N . 
Then 



\u\ 



\1 -LAj \P(v)(r'oP)(y,u)\ < \vr\y,u)\. 
An alternative formulation of the above in local coordinates is 



a I N _! 



dh dh 
dy dy 



a 1 1 



\u\ 



< \vi!(y,u) 



(75) 



(76) 



Proof. 

\vz(y,u)\ 2 = 



T I dh \ / dn \ i 2 



t , dh Ah ' 

v ■ lN -i + d^d^ 







T 
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Figure 4: Construction Of V\ 

„ / dh r dh\-" 2 i 

t! "Kd ( 



d 2 h 

dy 2 


T 









dn dn 

+ u 2 t> ( d y d y 
+ 2« I 1 





dhdh T \ 
dydy J 

dh T dh 
dy dy 
dh dh T 
dy dy 



J 



T 2 
3. -|- Clj\j 



-1/2 



d 2 /i T „ dndn T 
a— -a +ua- — — a 
dy^ dy dy 



a|/ 7V _ 1 + — ) a T [l + 2cm + /3u 2 ] + c4 



The proposition follows by applying Lemma [4.10 



□ 



We shall now remove references to local parametrisations of the boundary. 
If a point z G M. N is within a distance r of <9f2 then there is a unique nearest 
point of dfl to z. We shall define maps Ui, v-i which formalize the notion of 
a nearest point to dfl and Q respectively. Define V\ : (dfl) r — > dfl by 



-i 



v\ := pup 

This is well defined due to Lemma 14. 6L Define v 2 : fi r 



(77) 



n by 



V ; \^i(o;), io e p(<9ft x [0,r)). 
Note that p(dfl x [0,r)) n ft = because of the construction of p. 
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Figure 5: Contour sketch of d x for a typical region il 
Let iGll and define d x : Q r — > M+ by 

d^z) := d 9 (z, f 2 z) (79) 

— 2 

where we recall that the standard Riemannian metric d g : Vl — > M + is defined 
by 

d g (x,y) := inf{/(7) : 7(0) = x, 7(1) = y, 7 C f2, 7 is cts and piecewise C 1 }. 

(80) 

Construction 4.12. Let < 5 < r be fixed. Let 71 satas/y either of the 
following conditions, and define 72 accordingly. 

(i) Let 71 : [0, 1] — > (dQ)s be a C 1 curve such that 71 %. Q. Let 

T = M{t e [0,1] : 7l (t) (81) 

and /et 

T 1 = sup{tG [0,1] : 7l (t) (82) 
Define 72 : [0, 1] — > fry 

l2 {t) = [ Vlllitl tE ^ T ^ (83) 
1 7i(t), otherwise; 

(ii) Lei 7 i : [0, 1] — > Q be a C l curve. Then let 72 = 7 i. 
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Figure 6: Construction of the paths 71, 72, 73, and74 




Proof. We shall only prove the lemma for case (i), as case (ii) is trivial. For 
t G (T , Tx) let p G dfl be the nearest point of dQ to ji(t) and let (y, h(y)) 
represent the surface dfl in the local coordinate system based at p. Since 71 
is continuous, there exists an open interval / C (T ,Ti), containing t such 
that 7i(/) C t(G x (— r, r)). Let 73,74 be paths defined, for s G /, by 



Note that 7 2 |j = T74I7 and that 72,73,74 are C l on / because r is C 1 . By 
writing 73(3) = (a(s), ajv(s)) we see that |azv(£)| = d{pfi{t) , dVt) < 5. Thus 
by Proposition |4.11| , 



73 := T Si 



74 := ^73- 



(85) 



( 




< \mr'(ls(t))\ 



Integration yields the result. 



□ 



Corollary 4.14. 



(i) Let z\ gOj, where < 5 < r. Then 




(86) 
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(ii) Let x G fl and let z G -B(a;; 5), where < 5 < r. Then 

i i i \ 

d g (x, v 2 z) < 



1 1 

S r 



*7) 



(iii) Lei G f2 6e stzc/i t/iat d g (x,y) < 2r. Then 

2rd g (x, y) 



\y — x\ > 



2r + d g (x,y)' 

Proof, (i) Let z 2 G B(zi] 5 — d(zi, and let 71 be the straight line joining 
z\ and z%. Then 71 satisfies one of the two conditions in Construction |4.12 . 
Accordingly, 



14(^2) - 4(^0 



= \dg{X, IS 2 Z 2 ) - dg{X, V 2 Zl) 

< d g {y 2 z u v 2 z 2 ) 

< K12) 



\ r . 



1 - 



-1 



z 2 - Zx . 



(ii) Let 71 be the straight line joining x and z. Then 71 satisfies one of the 



conditions in Construction [4.12| . Hence, by Lemma |4. 13 



8 



5 



d g (x, v 2 z) ( 1 - -J < \1 - -J l( l2 ) < /(71) = 5. 

(iii) Let 71 be the straight line joining x and y. Then setting 5 = \\y — x\ < r, 
we see that 71 satisfies one of the conditions in Construction |4.12|. Hence, by 



Lemma 4.13 



1 



d g (x,y) <\y- x\. 



□ 



Rearranging this gives the desired inequality. 

Let ks be approximate identities, defined as follows: 

Definition 4.15. Let 5(0; 1) denote the unit ball in R N , and let k : 
-6(0; 1) — > R be smooth, non-negative and have unit integral. For 5 > 
define k s : B(0; 5) -> R by 



h{z) = 8~ N k (|) . 



59) 
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Let xefi and let (3 > K/r where 

K = K m , N>k : = sup ( I \D^k\) . (90) 

l<b'|<m-l \ J B(0;1) J 

Define f m>/3tX : Q -> R by 

Sm,p,x{y) ■= [ 1 - J / d x (z)k K/f3 (z - yjd^z 

= ( 1 -ir) I d ^ z + s/)**//^)^*- (91) 

\ /5r/ JB(0;K//3) 

Lemma 4.16. T7ie functions f m ,i3,x belong to the class & m p (see Defini- 
tion Ol). 



Proof. Suppose that 3/1,3/2 £ ^ and that |j| is a non- negative multi- index 
such that 1 < \j\ < m — 1. Then by Corollary |4.14j (i), 

, imsup PV^fe)-^7, g fa)l 

ie-n/i 1 3/2 — 3/1 1 

(! - f ) fB(0;K/p){ d *( z + 2/2) - + y 1 )}^'A; K//3 (2;)d Ar z| 



lim sup 



"1/2— vi 12/2-3/1 1 

< limsup / \df + yi)-d x (z + m)\ L_k\ 

yo^vi Jb(O;K/0) Z + Vo - Z + Vi ) \ Br T /PK ' 1 



B(0;K/f3) 



lB{o;K/p) \(z + y 2 ) - (z + yi)\ \ fir, 

limsup \ dx ^ z + y2) ~ d ^ z + Vi)\ L _ ^ 

V2-+yi \{z + y 2 ) - {z + Vl )\ V /3r, 



\D^k K/(3 (z)\d N z 



< \D j k K/0 (z)\d N z = K~^0 \D j k(z)\d N z < 0. 

JB(0;K//3) JB(0;1) 

Since f m ,p,x is C°° and satisfies the above inequality we see that | V f m ,/3,x(y)-u T \ < 
1 for all u € R N with |u| = 1, by taking j = 0, and hence that \V f m ,i3,x{y)\ < 
1. Also, letting j < i be any non-negative multi-index such that \j\ = \i\ — 1, 

|£>7^(2/)i<^ |ihl - □ 

Recall the definition of the Riemannian-type metrics d m ^ : Q 2 — > R + 

d m> p(x, y) : = sup{0(?/) - (f)(x) : G <f m>/3 } (92) 
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Lemma 4.17. Let f2 satisfy Condition [15]. Let x, y 6 Q and let (3 > K/r. 
Then 



d m ^(x,y) > d g (x,y) 1 



K\ 2K 



(93) 



Proof. If z e B(x;K/ (5) then by Corollary [4.14 



ii 



so 



d x (z) = dgix^z) < - J 



f ( K \ K 

frn,P,x( X ) = / d x{z) 1 - — k K/f3 (z - x)d N Z < —. 

JB(x;K/l3) \ (3r I p 



Similarly, if z G B(y; K/{3) then 



4(z) = d g (x, v^z) > d g (x, y) - d g (y, v^z) > d g {x, y) 



K r 



so 



K 



fm,(3,x(y) = / d A z ) 1 - IT} k K /f3(z ~ V)d Z 

JB{y;K/f3) \ (3r J 

( K\ K 



Thus using Lemma [4.16 



/ K\ 2K 

d m , p {x,y) > f m ,/3,x(y) ~ f m ,pA x ) > d g (x,y) I 1 - — J - — . 



□ 



The dependence of K upon k can be removed by taking the infimum of all 
values of K m)N>k) where k is smooth, non-negative, supported in the unit ball 
and has unit integral. For example, K 2) n — N 2 . 

Theorem 4.18. Let £1 satisfy Condition [O. F° r P > 4K/r we have 



1 - y j- T ) dg - dm ' p - dg ' 
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Proof. Let 7 C Q, 7(0) = x, 7(1) = y, 7 is piecewise C 1 and let be such 
that |V0| < 1. Then 

\<f>{y) - <f>{x)\ = |0(7(1)) - 0(7(0)) I 

;^(7(t))dt 

1 v<£.y (t)dt 



< / |V0||7'(t)|dt 

< i(7)- 

Taking the supremum over all <j) in S m ^^ we see that 



Let /3 > and let e := ^K/f3r < 1/2. For large distances d g (x 1 y), 

Lemma 4.17| is a useful result. We therefore consider two cases: 

(i) Suppose d g (x,y) < 2er/(l — e). Then we may use Corollary |4.14| (hi) 
to obtain 

d m ,p{x,y) > \y-x\ 



> 



> 



2rd g (x,y) 
2r + d g (x,y) 
2rd g (x,y) 
2r/(l-e) 



p,d g (x,y). 



(ii) For d g (x,y) > 2er/(l — e) we see that 



, , J K\ 2er 



e - e 



Hence 



2e 2 r 
2K 

T' 



, , J K\ 2K , , , . . 
dg{x, y) I 1 - — I - — > (1 - e)d g (x, y) 
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Using Lemma [4.17|, 



dm,fi(x, y)> ( 1 - \l — j d g (x, y). 



In both cases, we see that 



1 - d g( x >y) ^ d mA x >y) < d g (x,y). 



□ 



Note 4.19. Suppose that Q is an iV-dimensional locally Euclidean Rie- 
mannian manifold which posesses a locally injective isometric mapping into 
M. N (i.e. Q is a covering space of some non-simply connected open subset of 
M. N ). The notion of a Euclidean metric in such manifolds degenerates, and so 
heat kernel bounds involving such metrics are not useful. Heat kernel bounds 
involving the Riemannian metric can be found by adapting the result of this 
section and using the methods of Section 0. In fact, by using a covering 



space, it can be seen that Theorem |4.18| is valid where r is the greatest lower 
bound of the radii of curvature at all points of the boundary. □ 
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